Using electrowetting to control interface motion in patterned microchannels 
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We use mesoscale simulations to demonstrate the feasibility of a novel microfiuidic valve, which exploits 
Gibbs' pinning in microchannels patterned by posts or ridges, together with electrowetting. 



As fluidic devices become smaller there is a need for sim- 
ple and robust ways to control the motion of fluids moving 
through them. Several way to arrest capillary filling have been 
developed UHal- These are mainly based on using hydropho- 
bic patches |[fl-[3|] or narrow hydrophobic channels to arrest the 
flow iH3. In the first case filling is restarted by an electric 
field which decreases the contact angle, while in the second 
a driving pressure is applied. In this paper we use numerical 
simulations to motivate an alternative mechanism for creat- 
ing a microfiuidic valve. We show that placing obstacles at 
relevant positions along a microchannel, and addressing them 
by an electrowetting potential applied at the position of the 
interface 10], can allow capillary filling tU] to be halted and 
restarted reversibly. 

If a microchannel is brought into contact with a fluid, and 
the contact angle between the surface and the fluid, eq , is 
less than 90°, the fluid will move into the channel fl. This 
is the well-known process of capillary filling. However, if 
there are obstacles to the flow, such as micron-scale ridges or 
posts on the surface, the advancing front can pin for contact 
angles greater than a certain threshold value (0 tr < 90°), thus 
halting the filling. For a slowly moving interface the pinning 
can be understood in terms of the Gibbs' criterion [9] which 
states that an interface remains pinned to a corner as long as it 
makes an angle smaller than 6 cq with each of the dry surfaces 
meeting at the corner. As the interface is pulled along the 
channel by the capillary forces, it may not reach the required 
angle to wet the surface of the obstacles, and hence the filling 
is arrested. 

We consider first the geometry shown in Fig. QJ. Posts of 
square cross section DxD and spacing L reach across a chan- 
nel of height H. Flow is along the x direction and we consider 
a channel sufficiently wide that the effect of the side walls may 
be neglected. Fig. [Tj? depicts a pinned configuration for this 
geometry. Although the front is pulled along the top and bot- 
tom walls of the channel, for sufficiently large contact angles 
it does not move far enough to start wetting the dry face of the 
posts. For H/L > 1 the threshold contact angle below which 
flow is possible is dti ~ 55°. This value is independent of the 
channel height because the walls act independently IToll . As 
H/L decreases below unity, 9 tr increases to 90° as H / L — > 0. 
Hence a suitable choice of H/L allows #tr to be tuned as re- 
quired. 

Electrowetting [7], the application of an electric potential 
to a microchannel, is a robust method of changing the contact 
angles by tens of degrees. Hence one might expect that the ap- 
plication of an electrowetting potential can rather easily move 
the system between the pinned and the filling regimes, thus 
allowing control of the motion of the capillary front. In this 
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FIG. 1: (a) Channel geometry: the fluid flows along the x direction 
and periodic boundary conditions are used along the y direction, (b) 
Typical shape of the front when it is pinned at a row of posts. For 
clarity the top wall is omitted, (c) Position of the filling fluid versus 
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letter we present simulations of capillary filling in a patterned 
microchannel showing, in particular, how switching between 
contact angles above and below 6 tI can be used to halt and 
restart capillary filling. Before presenting our results we first 
outline the numerical approach we use. 

As we are considering micron length scales it is appropriate 
to describe the system using a mesoscale algorithm. We con- 
sider a binary fluid with components A and B, say, which in 
equilibrium is described by minimising the free energy func- 
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The subsequent position of the interface (recorded in the 
centre of channel) as a function of time is shown in Fig. QJ. 
The zero potential equilibrium contact angle of the channel is 
chosen to be 9 cq = 60°. This is greater than 8 t 



where n is the local total density of the A and B components 
(n = riA+ns), <fi is the order parameter tj> = ua—tib and c is 
the lattice velocity c = Sx/St, where Sx is the lattice spacing 
and St is the simulation time step. The first integral in Eq. (Q]i, 
taken over the total volume ft, controls the bulk properties of 
the system. The terms in <p give coexistence of phases with 
(\> = ±1. The energy cost of an interface between the two 
phases is modeled by the derivative term, with k related to the 
surface tension. The term in n controls the compressibility of 
the fluid. The integral over the solid-fluid interface dfl in Eq. 
(fl~|l accounts for the wetting properties of the solid surfaces 
111 111 . The surface field h can be used to tune the equilibrium 
contact angle. 

The hydrodynamics of the fluid is described by the Navier- 
Stokes equations for the density p and the velocity field v 
together with a convection-diffusion equation for the binary 
order parameter <f) 
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(1) the advancing front is pinned at the first row of posts, in the 
configuration reported in Fig.Q}>. At t = 2 ■ 10 5 St the contact 
angle of the channel is decreased to 8 V = 30°, corresponding 
to the application of an electrowetting field. Filling restarts 
immediately and the interface advances. It is accelerated (by 
the increased hydrophilic surface area) as it passes the second 
row of posts and then slowed as it reaches the trailing edges 
of the posts, but is not pinned 1 10]. We have checked that it 
also passes the third row of posts if 6 V is held at 30°. 

However, in the results presented in Fig. QJ-, the equilib- 
rium contact angle is switched back to 6 eq = 60° when the 
advancing front reaches the middle point between the second 
and the third posts. This leads to pinning, explicitly showing 
the possibility of stopping the imbibition of the channel at a 
prespecified position. 

As a more stringent test we also simulated a channel with 



d tP + \7-(pv)=0, (2) 
dt(pv) +v • V(pv) = -V • £ + ??V 2 v , (3) 
a t + V-(0v) =MVV (4) 

In Eq. ([3j 77 is the viscosity of the fluid and in Eq. (0]l M is a 
mobility coefficient. The pressure tensor P and the chemical 
potential p, which appear in Eqs. ([3]) and (|4J, are calculated 
from the free energy ([T). 

We have chosen to use a two-component binary fluid as 
a model system. This is because modeling capillary filling 
correctly using a liquid-gas diffuse interface model is com- 
putationally very demanding because of unphysical motion 
of the interface due to evaporation-condensation effects lfl2tl . 
However our results are equally applicable to a physical sys- 
tem where a liquid displaces a gas as the important physical 
parameters are the viscosities, not the densities, of the fluid 
components. Therefore we shall use the natural terminol- 
ogy 'liquid' and 'gas' for the displacing and displaced fluid 
from now on. The equations of motion are solved using a lat- 
tice Boltzmann algorithm [13]. Full details are given in Refs. 
lfl4l [l5ll . In particular lfl5ll shows that this approach gives ex- 
cellent agreement with analytic results for capillary filling in 
smooth microchannels. 

For our first simulations we investigate the microchan- 
nel geometry shown in Fig. Measuring distances in 

lattice units, the simulation box is of size (L x ; L y ; L z ) = 
(180; 60; 60) and we take periodic boundary conditions in the 
y-direction. The height of the channel H = 60, the separation 
of the posts L = 40, and the posts are of width D = 20 and lie 
between x = 20 and x = 40, x = 80 and x = 100, etc. At the 
beginning and end of the channel are liquid and gas reservoirs 
respectively, connected in order to equalise the pressure. The 
simulations commence with the interface at x — and with 
the elctrowetting potential switched off. 



six rows of posts switching between 6*e q = 60° and 6 V = 50°. 
Pinning and subsequent depinning when the electrowetting 
potential was activated was observed at each post. 

The broken line in Fig. [TJ: compares a similar simulation, 
but with the contact angle on the channel walls switched to 
9 V = 30° but that on the surfaces of the posts remaining at 
6 cq = 60°. The behaviour of the fluid front is essentially the 
same, except for a slight slowing due to the reduction in the 
capillary force. This demonstrates that the valve mechanism 
does not have a strong dependence on the details of the imple- 
mentation of the electrowetting potential. Moreover the ap- 
plied electric field does not need to be localised at a particular 
post (in our simulations we change the contact angle along all 
the channel); it is sufficient that it covers all the posts where 
filling is (re)-started by the applied potential. This avoids the 
challenge of miniaturising the electrodes. 

We next consider a second channel geometry, for compari- 
son, and to illustrate the robustness of the mechanism. This is 
shown in Fig.lJJj. The channel has height H = 50 and width 
W — 40 and side walls are now explicitly included. Only the 
bottom wall of the channel is patterned, by ridges of height 
and width 20 and spacing 40 which run all across its width. 
The aspect ratio of this channel, and the fact that only one wall 
is patterned, may make it more accessible experimentally. 

Running simulations to locate the threshold value of the 
equilibrium contact angle gives a value of 9 tr between 65° 
and 70°. A pinned configuration with 8 cq = 70° is reported 
in Fig. I2J7. (For a two dimensional geometry with no side 
walls 6* tr = 45°. Thus, as expected, the hydrophilic side walls 
significantly increase ease of filling.) 

Results of simulations for the grooved channel are shown in 
Fig. [2]:. The simulation procedure was analogous to that used 
for the channel with posts, but with the contact angle cycling 
from 6 cq = 70° at t = to 6 V = 55° at t = 2.2 ■ 10 5 St and 
back to 70° at t = 4.9 ■ 10 5 St. The pattern of filling is very 
similar; as before the interface is halted at the first and third 
posts (compare Figs.QJ- and Fig . EJ^) . 

If the electrowetting potential is applied to the top wall only, 
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FIG. 2: (a) Channel geometry: the fluid flows along the x direction, 
(b) Typical shape of the front when it is pinned at a ridge. For clarity 
the top and a side wall are omitted, (c) Position of the filling fluid 
versus time. In these simulations eq = 70°, 8 V = 55° and 8ti ~ 
65°. 

the front is still depinned but touches the subsequent post be- 



fore the bottom wall is completely wet. This leads to the 
formation of bubbles of air trapped between the posts. This 
could be avoided if necessary by choosing smaller ridges or 
more widely separated posts. An estimate of the geometry 
necessary to avoid bubble formation follows from observing 
that the advancing front makes an angle 6> cq with the dry face 
of the ridge as it depins. Therefore to avoid the front touch- 
ing the subsequent ridge without having completely wet the 
previous one, hence trapping a bubble, it is sufficient that the 
ratio of the height of the posts to the distance between them 
is less than tan(7r/2 — eq ). To check this we ran faster, two- 
dimensional simulations for the grooved geometry, similar to 
that reported in Pig. [2jz but without lateral walls and with the 
electrowetting potential applied just to the top wall. Keeping 
the equilibrium contact angle on the patterned substrate fixed 
at 60°, we switched the top wall between 60° and 25° repro- 
ducing a filling profile similar to that reported in Fig. [2]: and 
with no air entrainment. 

In this paper we have demonstrated that it is possible to 
use mesoscale simulations to study the behavior of an inter- 
face in a patterned microchannel. We have shown that surface 
patterning, together with electrowetting, can be used to con- 
trol how the interface advances along the microchannel. We 
do not expect the valve mechanism to depend qualitatively on 
the size of the channel as long as the channel is sufficiently 
small that gravitational effects can be ignored. The obstacles 
must be at least a few microns in size, so that the fluid can- 
not overcome the pinning due to thermal fluctuations, and the 
advancing fluid must be non-volatile or it will move past the 
obstacles by evaporating and then condensing. 

However there is a strong quantitative dependence of 8 tl - 
on the details of the patterning allowing this to be adjusted to 
lie between the available 8 cq and 6> v . Moreover, although we 
have considered the quasistatic case, the threshold angle will 
increase for an interface with finite velocity or for pressure 
driven flow. 

We hope that this paper will motivate experiments on pat- 
terned microchannels, and that simulations such as those pre- 
sented here will be useful in determining suitable channel ge- 
ometries for a given application and materials. Similar mod- 
elling approaches could also be used to predict when an inter- 
face depins from the end of a constriction in a microchannel 
or a micronozzle due to an applied pressure [6]. 
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